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Abstract 

The two-dimensional Peierls-Hubbard model is studied at half-filling within 
both Hartree-Fork and Kotliar-Ruckenstein slave-boson theory. The inter- 
play between two types of long-range order, bond-order wave (BOW) and 
antiferromagnetism (AFM), is analysed for two representative dimerization 
patterns, corresponding both to the same wavevector (vr, vr). For each pat- 
tern, the Peierls dimerization (and associated BOW) is weakened and finally 
suppressed with increasing Hubbard on-site interaction, and correspondingly 
AFM is gradually enhanced. In particular, a coexistence regime with both 
BOW and AFM order is established in the parameter space of electron-lattice 
and Hubbard interactions. 

PACS numbers: 71.45.Lr, 71.10.Fd, 75.30.Fr, 63.20.Kr 
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I. INTRODUCTION 



The Peierls instability towards spatially broken symmetry is an important phenomenon 
in low dimensional materials |l|]. The one-dimensional (ID) case has been widely discussed 
in the context of polyacetylene (CH)^: based on the Su-Schrieffer-Heeger (SSH) model [0,0, 
where lattice displacements couple to electron hopping. For a half-filled band arbitrary 
small electron-lattice (e-1) coupling will induce a lattice dimerization (disregarding quantum 
lattice fiuctuations), which is associated with a periodic modulation of the bond hopping, 
that is, a so called on-bond charge- density wave or bond-order wave (BOW) [Q. It has 
been established that the Hubbard on-site Coulomb electron-electron (e-e) interaction U 
will enhance the bond alternation initially for small values and finally suppress it at large 
values of U 

In two dimensions few theoretical investigations exist [[7|-|Tl[|, some of which connect 
the physics of Peierls systems to that of the high-T^ copper oxides Moreover, these 

investigations may be of direct relevance to those quasi-two dimensional (2D) materials which 
show a Peierls instability such as transition-metal oxide bronzes like AMoeOiy (A=Na, K, Tl) 
g and organic conductors like (BEDT-TTF)2MHg(SCN)4 (M=K, Rb, Tl) As an 



effective minimal model, in this context, the 2D version of the SSH model was investigated 
PJ5|jTI|]. With only nearest-neighbour (n.n.) hopping on a square lattice, the electronic 
Fermi surface is perfectly nested at half-filling with nesting vector Q = (tt, tt). Two possible 
alternation patterns for the lattice distortion and the concurrent bond hopping comply 
with this Q, as illustrated in Fig. |l[ Whereas for pattern (a) the dimerization is in both 
directions, it is only in one direction for pattern (b) [§. Similar to the ID case, already for 
an infinitesimal e-1 coupling, the 2D SSH model goes through a Peierls instability into one 
of the dimerized states of Fig. [l| pJ]) . 

When a Hubbard on-site Coulomb interaction U is included — the model is then the 
so called Peierls-Hubbard model, results differ from the ID case. Numerical calculations 
on a small 2D lattice MM indicated that the Peierls instability will be suppressed by the 



2 



Hubbard U [|T^. An intuitive explanation is that the on-site Coulomb interaction favors a 
spin-density wave (SDW) long-range order, that is, antiferromagnetism (AFM), while the 
dimerization associated with BOW harmonizes with a spin-singlet formation between those 
two n.n. spins which are connected by a strong bond. As we know, due to the same nesting 
effect, the pure 2D half-filled Hubbard model (without consideration of a Peierls instability) 
has been shown to exhibit AFM long-range order for any U > 0. This is in stark contrast 
to the corresponding ID case where no true long-range AFM order is available and the 
correlated state rather corresponds to a resonant valence bond state with strong weight 
from n.n. singlets [jl5|. Consequently, one may envisage, for finite e-1 coupling (denoted 
as 7], see below) and e-e on-site interaction U, a competition between BOW and AFM as 
the underlying physics in the 2D half-filled Peierls- Hubbard model. In the large U limit, 
Zhang and Prelovsek have studied the corresponding spin-Peierls (SP) instability and found 
that the SP state, competing with AFM, does not appear unless the spin-lattice coupling 
(analogous to r] here) exceeds a threshold |T^. 

The details of the competition between the two ordered states were studied only for the 
above limiting case and the situation is not clear for general values of t] and U. In particular, 
a basic problem has to be solved: does BOW disappear once the AFM order sets in, or is a 
coexistence of the two orders possible? It was previously argued by Mazumdar within a real 
space approach that the appearance of the AFM should coincide with the disappearance of 
the BOW ip, which was, however, not verified. To clarify this issue — which is the topic 
of this paper, one needs to explicitly calculate the two order parameters for the BOW and 
AFM with varying rj and/or U. 
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FIGURES 




(a) (b) 

FIG. 1. The lattice distortion patterns (a) and (b). In the figure a thick sohd Une corresponds 
to a strong bond with hopping integral t(l + 6), a dashed line corresponds to a weak bond with 
hopping integral t{l — 5), and a thin solid line corresponds to a normal bond with hopping integral 
t. Both patterns correspond to phonons with wave vector (vr,7r). The dimerization is along two 
axes for pattern (a), while only along the x axis for pattern (b). 

In this paper we make use of both the Hartree-Fock (HF) and Kotliar-Ruckenstein slave- 
boson (SB) approach []17| to treat the Hubbard interaction. The HF results are usually 
valid at weak-coupling, and they can be used as a basis for further elaborate studies. In the 
context of investigations on density wave instabilities, the HF theory may give reasonable 
results even in one dimension [jl8|, where one would expect it to be worst because of strong 
fluctuations. In dimensions higher than one, as considered here, qualitatively correct results 
are expected from the HF theory. In order to extend the controlled weak-coupling results 
to intermediate values of U, we evaluate BOW and AFM within a slave boson mean-field 
approach, which is considered to be appropriate to interpolate from weak to strong coupling 

The paper is organized as follows. In the next section the model Hamiltonian is intro- 
duced, and then the HF and SB approaches are formulated. The self- consistent equations 
for the order parameters are derived in both theories. In Sec. Ill numerical results are 
presented. The main results are shown in Fig. ^, where the coexistence of BOW and AFM 



is found to be possible for each of the two patterns, in contrast to Mazumdar's argument. A 
complete comparison is made between the results derived from both approaches. Conclusive 
remarks are given in Sec. IV. An appendix completes the SB formulation. 

II. FORMULATION 

We begin with the 2D half-filled Peierls-Hubbard model 

H = H, + Hu + Hk (1) 

with 

Ht = -tY.\^ + "K, - <+l,i)](4j,aCi+lj,a + h.c.) 

-t ^[1 + a{ulj - Mf_j-+i)](4j-^Qj+i,^ + h.c.) , 
^ hi 

where c|j ^(cjj^o-) is the creation (annihilation) operator for an electron at site («, j) with spin 
a {i denotes x coordinate and j denotes y coordinate), njj^o- is defined as Ui^j^a- = cl,j,crCj,i,o-) 
u^Jj is the displacement component of site (2,j) in x/y direction, t is the n.n. hopping 
parameter, and a is the electron-lattice coupling constant. Hjj is the Hubbard on-site 
interaction with the repulsion strength U . The last term Hk is the lattice elastic potential 
energy, with K the elastic constant. The phonons are treated in adiabatic approximation. 

For an analytical treatment on an infinite lattice, we have to work with a definite distorted 
lattice, rather than allowing the distortions to arise arbitrarily. In this paper we constrain 
the discussion to the lattice distortions within the two commonly used dimerization patterns 
shown in Fig. |l|. Only these patterns correspond to the nesting vector Q = (vr, vr) and they 
realize an unconditional Peierls instability that occurs for a — > and U = 0. Explicitly they 
are written as 

<, - = i-iy-'n, - = i-iy^^u 

5 



for pattern (a) and 

for pattern (b). For convenience, two dimensionless parameters are defined: the dimerization 
amplitude 5 — au and the electron-lattice coupling constant rj — aH/K. Throughout the 
paper the hopping integral t is taken as the energy unit. 

In the following we will construct the analytical formulas based on the HF and SB 
approaches, respectively, and leave the numerical calculations to the next section. 

A. Hartree-Fock theory 

In our model, the on-site charge density wave is not favored and the total electron number 
on each site is uniform and equal to one at half-filling. Then the expectation value of the 
electron density with a given spin may be simply assumed as {rii^j^a) = |[1 + cr^—iy^^m] 
when the AFM order is taken into account, where m represents the staggered magnetization. 
In HF approximation (equivalent to Hartree here) the local Hubbard term may be decoupled 

as t/Eij ^ij.T^ijU ^ ^Sij(riij,T(^ij,i) + (wi,j,T)^i,ja-(^ij,T)(^iju))- Then the Hamiltonian 
becomes quadratic and may be easily diagonalized in momentum space. Under consideration 
of a bipartite lattice the final electronic spectra are derived as follows for pattern (a) and 
(b), respectively, 

£j„ = ±^JUV/A + 4[(cos fca; + cosy 2 + 52(sin^^ + sm kyf] , (2) 
= ±^Ju'^m? /4 + 4 [(cos K + cos kyf + 5'^ sin^ A;^] . (3) 

Each branch above (— or -|-) is two- fold (spin) degenerate. The wave vector \i — [k^, ky) 
is restricted to the reduced Brillouin zone: — tt < kx ± ky < tt. With inclusion of constant 
terms the ground state energy is 

E,^2J2 + NU{1 + m^)/4 + El,u , 

k 
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where i/ = a, b represent pattern (a) and (b), respectively, N is the total number of lattice 
sites and El^^, denote the lattice elastic energies for both patterns: E^^a = '^E^ b = NS^/rj. 

The self-consistent equations for dimerization 6 and magnetization m are found by min- 
imization of the ground state energy. They read: dEy/d5 = and dEy/dm = 0. The latter 
results in (except for a trivial solution m = 0): 



and the former leads to 



, 4?7 ^ (sinfc.^ + sin ky)"^ . 



for pattern (a) and 



N k Ka 



8r] ^ sin^ k^ . 

k l^k.fel 



for pattern (b). 

B. Slave-Boson theory 

In the spirit of the Kotliar-Ruckenstein SB approach [O, four auxiliary bosons 



ejj, pl^j^^{a =t,i), c^lj- are introduced to label the four different states for an arbitrary 
site which can be empty, singly occupied by an electron with spin up or down, or 

doubly occupied. The unphysical states in the enlarged Hilbert space are ehminated by 
imposing two sets of local constraints: 

^ij^iJ + J2plj,aPi,j,a + 4,jdij = 1 (completeness) , (7) 

and 

pl,j,crPi,j,a- + d\ jdi,j = clj^^Cij^o- (correctness of fermion occupancy for a given spin) . (8) 

For a bipartite lattice, we introduce a set of bosons, with separate Lagrange multipliers for 
each sublattice. At the mean-field level, the bosons are replaced by c-numbers and assumed 



to be site-independent on each sublattice. At the same time, the constraints above are 
softened to be satisfied only on the average on each sublattice. This treatment is equivalent 
to making a saddle-point approximation in the path-integral formulation. For concreteness, 
we introduce the following parametrization for sublattice A (and similar parameters are 
defined for sublattice B): e^, PAa, dA as average values of the boson operators e-|j- , p^i],„, d[^], 
respectively, and \a, as Lagrange multipliers associated with the constraints (0), (||), 
respectively. Then the Hamiltonian may be recast into the following form [we choose 
pattern (a) as an example]: 

H = -t{l + 6) J2 ^<^i4,j,abi+ld,a + + + h.c.) 

Ri,j€A,a- 

-t{l -6) zM,j,aai+l,3,a + ^a^j+i,, + h.c.) + NU{d\ + 4)/2 + EL^a 

-Ab E (4 + EpL + 4 - 1) - E Ab(pL + 4 - hl^M,,^) , (9) 

where a^(a) and 6^(6) are the electron creation (annihilation) operators for sublattice A and 
B, respectively. The hopping renormalization factor z^r ensures the correct result in the limit 
of vanishing U and takes the form Za = {z^){z^) with 

{z'^)= , -LPL.+PL.d, L = A,B. 

The Hamiltonian (P) may be diagonalized in momentum space and the energy bands read: 



4a,a = + A^)/2 ± V(A^ - A^)V4 + 4^2 [(cos fc, + COS A;,)2 + 6^{smK + sin A;,)^] . (10) 
Similarly, the energy bands for pattern (b) are given by: 



4a,b = i^A + A^)/2 ± V(A^ - AS)V4 + 44[(cos A;. + cos A;,)^ + sm' K] . (11) 

At half-filling and zero temperature only the two lower (— ) bands are occupied (the constant 
Aa + Ab is independent of a as will be seen later). Then the ground state energy is expressed 
as {u = a, b) 

8 



^^ = EW + ^0 + i?L,. (12) 
ko- 

with the constant Eo = {N/2)[U{d\ + 4) - XA{e\ + EaPla + 4 - 1) - Ai?(e| + EaPL + 

4 - 1) - A^(pL + 4) - A^(pL + 4)] • 

The self-consistent equations are obtained from the requirement that the ground state 

energy is stationary with respect to the parameters: eA(B), PAiB)a, dA(B), Xa{b), ^A{B)y ^■ 
Except for the equation corresponding to 6, they all have the general form: J2ka ^^k(Ti//^^ + 
dEo/dX = 0, where X represents one of the parameters. Analyzing these equations and 
applying the constraint Y.aP\{B)u + ^^(b) = Eo-(a-(&)Jj,o-'^(^)j,j,cr) = 1 at half-filling, one 
may find the solution satisfying the following relations: = = c^a = '^b(= d), \a = 
^b{= A), pAa = PBa, = ^B Ect ^a — Ecr ^B = ^- Consequently, the number of free 
parameters is substantially reduced. The final compact self-consistent equations are listed 

t t 

in the Appendix with several re-defined independent parameters: d, A, Xab = ^a ~ ■^b = 
^B ^ A;^, m = — = — where m denotes the same staggered magnetization 
as in HF theory. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

We now focus on the numerical results obtained from the self-consistent equations. First, 
it is necessary to analytically analyze the mean-field equations more thoroughly. After a 
replacement of the momentum summation by integration, i.e., X^k ~^ ^ I dfcid/c2 
{kx,y = ±^1 + ^2), we examine Eq. @ in HF theory. It is easily seen that the right-hand side 
(rhs) of Eq. (Q) assumes a different analytical behavior for each of the two patterns due to 
their different respective spectra. For pattern (a) the rhs is divergent at m = 0, irrespective 
of the value of 6, and decreases monotonically with increasing m until it reaches a value less 
than 1 at m = 1. This implies for pattern (a) that Eq. (^) can be always solved with a 
nonzero solution of m as long as f/ > 0. For pattern (b) it is not always the case because the 
rhs is finite at m = for any 6 > 0. Once this finite value is less than 1, one has to adopt 
the trivial m = solution. The same conclusion may be obtained by a similar analysis of 
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the corresponding equations in the SB evaluation [see Eqs. ( |X3| ) and (^) in the Appendix]. 
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FIG. 2. The optimal values for the dimerization 5 and AFM m as functions of [/ at ?? = 0.5. 
In each panel, the dashed lines are for pattern (a) and the solid ones are for pattern (b); the thick 
lines are the SB results and the thin ones are the HF results. 

The numerical solutions for 5 and m are displayed as functions of U in Fig. ^ for both 
patterns and in both approaches, whereby the e-1 coupling rj is fixed at 0.5. Globally, it is 
seen that for both patterns 6 tends to decrease and m to increase with growing U . This 
supports the notion that the on-site interaction tends to favor the AFM order and to suppress 
the Peierls dimerization. Let us go to the details in the following. 

In HF (see all the thin lines), it is found that for pattern (a) m becomes finite (although 
small at small U) and simultaneously 5 begins to decrease once [/ > 0; while for pattern (b) 
m stays zero for small U up to t/ > Um — 1.34 where it becomes finite and correspondingly, 
5 first keeps its t/ = value and then starts to decrease for U > Um- The dimerization 
disappears at the same critical value Uc — 2 for both patterns and it approaches zero 
smoothly and quickly when U is close to Uc- Comparing the results for m vs. U between the 
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two patterns, it is clear that pattern (a) is more favorable to the formation of AFM order 
than pattern (b). We will come to this point later. 

Most of the above qualitative results are also found in SB approach (see all the thick 
lines). On the other hand, the difference to the HF results is clear as well, which we 
want to emphasize here. A distinct quantitative difference is that for each of the patterns 
the AFM order derived from SB theory is (much) weaker than that from HF theory, and 
correspondingly the dimerization decreases to zero over a larger range of U. The critical value 
for the disappearance of 6 now becomes Uc — 2.38, and the necessary Hubbard interaction 
to induce finite m in the case of pattern (b) is Um — 1.56. Both values are larger than the 
corresponding ones from the HF theory, which is understandable. As is well known, the 
HF theory usually overestimates the tendency towards AFM order. The SB theory, as an 
improved approach to fluctuation contributions, should lead to a slower formation of AFM, 
which complies with our findings. 

A further important difference between the SB and HF results is observed in the region 
U < Um for pattern (b), where the AFM order has not yet formed. From the solid lines in 
the upper panel of Fig. 0, it is seen that the dimerization 6 keeps a constant value in HF 
theory, while it decreases slowly with increase of U in SB theory. This disagreement may be 
understood as follows. In HF theory, see e.g., Eq. (|]), the Hubbard U becomes irrelevant 
once the order parameter m is zero: the value for 6 will be the same as that without U. 
In the SB approach, however, the Hubbard U is relevant even at m = by affecting the 
probability of double occupancy (P. As U increases, the double occupancy is disfavored, i.e., 
the quantity (f decreases. Correspondingly the effective hopping tz^- decreases too (cf. the 
expression for W in the Appendix), which may be understood equivalently as an increase 
of the elastic constant i^' or a reduction of the e-1 coupling rj. This signifies a decreasing 
dimerization. In this point the HF theory fails to catch the correct physics by assuming the 
probability of double occupancy as a constant 1/4 which is correct only for U ^ 0. 
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FIG. 3. The ground state energies (per site) as functions of f7 at 77 = 0.5, corresponding to 
Fig. |2| (with the same hue labels). Each line stops at its own critical point Uc- The original values 
of these energies are shown in the inset. All energies are in units of t. 

Furthermore, we can check the stabihty between both patterns by comparison of the 
ground state energies calculated in all cases, which are shown in Fig. ^. For each pattern 
the SB approach gives a lower energy than the HF theory in the whole range of U. Also, 
it is seen within each approach that pattern (b) has a lower energy than pattern (a), which 
signifies that pattern (b) is more stable. 

It is worthwhile to point out that the BOW is always associated with the finite dimer- 
ization. In order to see this, we have calculated the BOW which is characterized by 
a modulation of the hopping amplitude. Explicitly, we define the expectation values 

= {alj^^bi+ij^o- + h.c), h2 = (fej j^o-O-i+ij.o- + h.c.) {a is irrelevant) for the alternating 
bond hoppings along the x axis, and similar values h^, h\ for the hoppings along the y axis. 
By symmetry, we have — ^12 pattern (a) and h\ = h\ for pattern (b). All the 
quantities are calculated in both theories and plotted as functions of U in Fig. 0. The BOW 



is exhibited by the inequality between h\ and /if for each pattern ||20|. It is clear for each 
approach that such an inequality is present within < f/ < Uc, the same region where the 
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dimerization is finite. 
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FIG. 4. The expectation values /if 2) /if 2 ^ functions of U for both patterns. In each panel, 
the solid lines show the quantity /if and the dot-dashed lines give /if; the thick lines are the SB 
results and the thin ones are the HF results. For pattern (a), /if 2 = /if 2- -'^^'^ pattern (b), h\ = /12, 
which are plotted by the dashed lines in the lower panel. 

Numerically, Tang and Hirsch [§] studied the same model and calculated the energy gain 
from dimerization as a function of U for the pattern (b) shown here. By studying how the 
energy gain changes with f/, they found originally that the Hubbard U has little effect on the 
dimerization until it is large enough to suppress it, and later corrected that the dimerization 
is disfavored as soon as U is present. The finite-size effect was cautioned by the authors 
themselves. Their principal result, i.e., the Hubbard U is unfavorable to dimerization, is 
consistent with ours, especially with the SB results for pattern (b). Although it seems that 
the suppression of dimerization by U is faster in our results than what they displayed, no 
direct comparison is available because they calculated neither the order parameters nor the 
ground state energies. Obviously, further numerical calculations on large size systems are 
necessary for better comparison. 
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FIG. 5. Stable phases for BOW, AFM, and the coexisting state (C) in the parameter space 
{U,rj) for pattern (a) [upper panel] and pattern (b) [lower panel]. The thick lines show the SB 
results and the thin ones show the HF results. 

The main contribution in our work is that the order parameters as functions of U are 
explicitly obtained so that the interplay between BOW and AFM becomes transparent. The 
problem proposed in the Introduction is then naturally answered. It is clearly seen in Fig. |^ 
that the BOW and AFM may coexist for both patterns. For pattern (a) the coexistence 
(i.e., 6 > 0, m > 0) appears as long as f/ > and pure AFM order exists for U > Uc- 
For pattern (b) the coexistence is limited to the region Um < U < Uc- These results are 
not favorable to the argument by Mazumdar that the AFM emerges with the disappearance 
of the BOW 0. In fact, the valence bond approach adopted by Mazumdar in real space 
is appealing. It states that, in order to implement a symmetry-broken state (e.g., BOW), 
"extreme configurations" with shortest repeat units have to be identified. For 2D systems, 
he chose pattern (a) as the realization of the Peierls state and argued that the extreme 
configuration for BOW is a combination of zigzag chains and the n.n. sites within each 
chain are doubly occupied and unoccupied, respectively. However, the pattern selected in 
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his work is not the pattern with the lowest energy and furthermore, the considered extreme 
configuration is actually disfavored, even for weak Coulomb on-site interaction, as is verified 
in the exact U approach of this paper. It implies that spin fiuctuations are more 
pronounced than charge fiuctuations in the BOW state of the half-filled system. 

Eventually, we determine the coexistence regions for different rj in both theories. The 
results are shown in Fig. |^, where different phases are indicated in the parameter plane 
{U, rf). For pattern (a), only two phases exist, either a state with coexisting BOW and AFM 
or a pure AFM state. However for pattern (b) pure BOW and AFM states exist, which are 
separated by a coexisting state — the region between two thick (SB) or thin (HF) lines in 
Fig. 1^. As for the methods, globally speaking, the SB approach pushes the AFM order to 
the higher U regime than the HF theory. 

Finally, we come back to the difference between the results for the two patterns. As 
explained above, pattern (a) is more favorable to the development of AFM than pattern (b). 
It may be roughly understood from their different dimerization structures. As seen from 
Fig. for pattern (b), each site is connected to only one n.n. site by a strong bond when 
the BOW (or dimerization) forms. Thus a spin singlet on this strong bond is apt to prevail 
in presence of f/, which will resist the AFM. On the other hand, for pattern (a), each site 
connects two n.n. sites with strong bonds. This, on the contrary, makes the construction of 
spin singlets on these strong bonds difficult and the AFM is easier to develop. 

IV. CONCLUSION 

We have investigated the Peierls-Hubbard model in two dimensions at half-filling within 
both HF and SB approach. Two dimerization patterns, corresponding to the same wavevec- 
tor (7r,7r), are considered and the interplay between two long-range order states, BOW and 
AFM, is addressed. For each pattern, it is found that the Peierls dimerization (and asso- 
ciated BOW) is weakened by the on-site interaction U as soon as U is present and finally 
suppressed at a critical U = Uc- Correspondingly, the AFM is favored by U. Whereas for 
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pattern (a), see Fig. AFM is induced once [/ > 0, it is not activated until U = Um for 
pattern (b). For both patterns, the coexistence of BOW and AFM is possible. SB and 
HF evaluations lead mostly to the same qualitative results and quantitatively the former 
approach results in larger values of Uc and Um- Whereas the HF evaluation provides us with 
the exact weak coupling {U 0) result, the SB approach extends the findings to intermedi- 
ate coupling, and corrects charge and spin fluctuation contributions beyond HF. Especially, 
the reduction of charge fluctuations by U decreases the dimerization 6 consistently in the 
region U < Um for pattern (b). 
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APPENDIX: SELF-CONSISTENT EQUATIONS IN THE SB THEORY 

In this appendix we implement the formulation of the SB theory. With respect to the 
parameters d, A, Xab, the self-consistent equations are derived as follows for pattern 

(a): 

(Al) 
(A2) 
(A3) 
(A4) 
(A5) 

where 



d 


= -2C1/2/A , 










-A 


= [7/2 + 2(^2/^/(14 


-m)/2 






m)/2-rf2)/2 


Xab 


= 4(C3/v/(l-m)/2 


-d^- 




\l + m)/2 




m 


= —XabIi , 










1 


= ^vWh , 
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h = 




h = 


1 


l3 = 




w = 


zl-- 



1 

+ 4iy[(cos /ca; + cos/cj/)^ + (52(sin /c^; + sin/cj^^"^ 
(cos + cos /ej^)^ + 5^(sin k^ + sin fc^^)^ 
+ 4Vr[(cos /^a; + cos kyY + 52(sin fc^ + sin kyY\ 
(sin fc^ + sin kyf' 
\\bI^ + 4W^[(cos k^ + cos A;j^)2 + ^^(sin k^ + sin kyY] 
lQd\^{l + m)/2 - d2 + ^(1 _ m)/2 - d?Y 

(1 — iv?Y ' 

64d3(y^(l + rn)/2-(P + ^(1 - m)/2 - ci^)^ 
(1 — m?Y 

'l-m)/2-<P[{l - m)/2 + d^] + (P J {I + m) /2 - <P 

'■ V m ^ —m \ , 

1 — m ' 



64d\,J{l + m)/2-d^ + ^'(l -m)/2-rf2)3(l - rf2 + ^'(l + m)/2 - d^^{l - m)/2 - d^) 

(1 — m2)2(l — m) 

Ud\^{l+m)/2-d'^ + y^(l -m)/2-rf2)3(l - rf2 + y^(l + m)/2 - rf2^(l - m)/2 - rf2) 



(1 — m2)2(l + rn) 

The equations for pattern (b) are the same as those above except that the expression (sin kx + 
sinfcj^)2 in /i,2,3 is substituted by sin^ fc^., and the Eq. ( |A5[ ) is changed into 

1 = ST]Wh ■ (A6) 

Correspondingly, the ground state energies may be written in the simple form: 



Ea = -2Y,\J>^ab/^ + 4iy[(cosA;a. + cosA;j^)2 + d'^^sink^ + sinky)'^]+NUd'^-Nm\AB/'^+EL,a 

k 

for pattern (a) and 



Efe = -2 ^ ^\\b/^ + 4iy [(cos fc^ + cos kyY + 52 sin2 k^] + NUd^ - NmXAB/2 + EL,b 

k 

for pattern (b). 

We point out that for both patterns there always exists a trivial solution with m = 
[consider Eqs. ( |A3| ) and ( |A4| ) and note that C2 = C3 at m = 0]. Moreover, we checked that 
for 5 = the results presented by Fresard et al. (e.g., staggered magnetization, ground state 
energy) for the pure Hubbard model ||l9l are reproduced. 
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